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Abstract 

Superkamiokande atmospheric neutrino data for 414 days are analyzed in 
the framework of three flavor oscillations with mass hierarchy. It is shown that the 
best fit point is very close to the pure maximal «-> v T case and Am 2 ~ 7 x 10~ 3 
eV 2 . The allowed region at 90 %CL is given and the implications to the long 
baseline experiments are briefly discussed. It is also shown that the threefold 
maximal mixing model fits to the data of atmospheric neutrinos and the reactor 
experiments for Am 2 ~ 1 x 10~ 3 eV 2 . 

I. Introduction 

Recent atmospheric neutrino data of the Superkamiokande experiment 

II, 12,17 proving very strong evidence for neutrino oscillations. The original analysis 
by the Superkamiokande group is based on the two flavor framework. Despite the 
constraints from the reactor experiments 27 ' 1 ' 2 , there can be small mixing between 

and v e in the atmospheric neutrinos. To estimate the oscillation probability 
of v e <-> v x or v e «-> v x in the long baseline experiments 23 ' 21,16,25,26 , therefore, it is 
important to know the allowed region of the mixing parameters at a certain con- 
fidence level. In this talk I present a three flavor analysis of the Superkamiokande 
atmospheric neutrino data for 414 days 11 . Some of the contents of this presenta- 
tion overlaps with Lisi's talk 18,7 . 
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2. Three flavor analysis of the Superkamiokande atmospheric neutrino 
data 



To evaluate the number of events, we integrate numerically the Schrodinger 
equation 
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is the MNS mixing matrix 19 , E is the neutrino energy, A = ^/2GpN e (x) stands 
for the matter effect 20 ' 28 in the Earth, AEij = Am|/2£, = mj — rrij. We 

assume without loss of generality that Am 2 x has the smallest absolute value of 
the mass squared difference. To account for the solar neutrino deficit 3 we assume 
that Am^ is of order 10" 5 eV 2 or KT 10 - l(r 9 eV 2 . On the other hand, the zenith 
angle dependence in the atmospheric neutrino anomaly requires that the larger 
mass squared difference be of order 10~ 3 — 10~ 2 eV 2 . So the pattern of the mass 
squared differences has hierarchy and can be classified into two cases which are 
depicted in Fig. 1. 
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Fig. 1 The hierarchical neutrino mass squared differences. The scenarios 



(a) and (b) are related by exchanging m 2 <-> — m 2 
in vacuum but physically inequivalent in matter. 



They are equivalent 
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The eigenvalues of the Schrodinger equation (|1|) are given by 4 




where 
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Under our assumption that Am 2 < 10~ 5 eV 2 , in the contained atmospheric neu- 
trino events (0.1 GeV < E < 100 GeV) we have \AE 2 \L\ <C 1 for downward 
going events (L is the neutrino path length and L< 500 km), and |A£?3i|, 
|A_E 32 |, A ^> |A£?2i| for upward going events. Thus the eigenvalues in (g) 
with Am^ ~ O(10~ 5 eV 2 ) or (9(10~ 10 eV 2 ) are almost the same as those with 
Am^ = to good precision, i.e., the matter effects are much more important than 
the smaller mass squared difference Am 2l < 10~ 5 e V 2 in the atmospheric neutrino 
problem. Hence we can safely ignore Am 21 in the following discussions. In the 
limit Am 21 = the mass matrix on the right hand side of (P becomes 
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where D = diag (e lS , 1, lj, we have used the standard parametrization 24 of the 
MNS mixing matrix 19 
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are the Gell-Mann matrices. As we see in (Q), 6*12 drops out and so does 5, as 5 
appears only in the overall phase of the oscillation amplitude A{y a — > z/g). Putting 
m 2 = Am§ l5 therefore, we are left with the three parameters (m 2 , 9\ 3 , #23)- 
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The way to obtain the numbers of events is exactly the same as in Foot et 
al. 8 , and we refer to Foot et al. 8 for details. In Foot et al. 8 two quantities have 
been introduced to perform a x 2 analysis. One is the double ratio 15,10 



R = 



(Wlc 

(iVi\y| no 



where the quantities N ejfM are the numbers of e-like and //-like events. The nu- 
merator denotes numbers with oscillation probability obtained by (|l|), while the 
denominator the numbers expected with oscillations switched off. The other one 
is the quantity on up-down flux asymmetries for a-like (a=e,/z) events and is 
defined by 

V _ (A^ViV+^losc 
* rv 



— (iV -o. 2 /iV+o. 2 )|„' 

where N~ v denotes the number of a-like events produced in the detector with 
zenith angle cos 6 < — r], while N+ v denotes the analogous quantity for cos > rj, 
where r\ is defined to be positive. Superkamiokande divides the (—1, +1) interval 
in cos into five equal bins, so we choose rj = 0.2 in order to use all the data in 
the other four bins. Thus x 2 f° r atmospheric neutrinos is defined by 
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where the sum is over the sub-GeV and multi-GeV cases, the measured Su- 
perkamiokande values and errors are denoted by the superscript "SK" and the 
theoretical predictions for the quantities are labeled by "th". In Foot et al. 8 a x 2 
analysis has been performed using the quantities R and Y's, or using Y's only. 
Throughout this paper we use the quantities R and Y's to get narrower allowed 
regions for the parameters. We have to incorporate also the results of the reactor 
experiments. We define the following x 2 '- 

CHOOZ / \ 2 Bugey / _ \ 2 Krasnoyarsk / _ \ 2 

j=X,12 \ UX J / i=l,60 \ ux ? / j=l,8 \ UJj 3 / 

where Xj are experimental values and are the corresponding theoretical pre- 
dictions, and the sum is over 12, 60, 8 energy bins of data of CHOOZ 2 , Bugey 1 
and Krasnoyarsk 27 , respectively. There are 6 atmospheric and 80 reactor pieces of 
data in x 2 = Xatm + Xreactor an d 3 adjustable parameters, m 2 , 6*13 and #23, leaving 
83 degrees of freedom. 

The best fit is obtained for (m 2 , tan 2 6> 13 , tan 2 623, x 2 ) = (7 x 10~ 3 eV 2 , 
1.1 x 10" 2 , 0.93, 71.8) for m 2 > and (-7 x 10" 3 eV 2 , 1.5 x 10~ 2 , 1.1, 71.4) 
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for m 2 < 0, respectively. xLn = (Xatm) mi „ + (xLctor) mi „ = 4.2 + 67.2 = 71.4 
indicates that a fit to data is good for 83 degrees of freedom at the best fit point. 
The allowed region for m 2 with 13 , 9 2 3 unconstrained is given in Fig. 2, where 
A X 2 = Xatm + xLctor - (Xatm + xLctor) mi „ < 3.5, 6.3, 11.5 corresponds to la, 
90 % CL and 99 % CL, respectively. The allowed region for \m 2 \ at 99% CL is 
4 x 10" 4 eV 2 < \m 2 \ < 1.7 x 1(T 2 eV 2 . 
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Fig. 2 

Value Of Ax 2 = Xatm+X^actor- 
(Xatm Xreactor ) min Xatm 

xLctor-71-4. The solid, dash- 
dotted, dotted lines represent 
the scenarios (a), (b) and the 
two flavor case with maximal 
v u, v T mixing, respectively. 



Using the same parametrization as that in Fogli et al. 5 , the results for the 
allowed region of the mixing angles (9 13 , 6 2 z) are given for various values of m 2 in 
Figs. 3 and 4. The results for m 2 > (Fig. 3) and m 2 < (Fig. 4) are almost 
the same. Notice that #i 3 can be large in the region of m 2 where the CHOOZ 
data 2 give no constraint on 9 13 . 
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Fig. 3 The solid, dashed, dotted lines represent the allowed region at 68 
% CL, 90 % CL, 99 % CL, respectively for scenario (a) in Fig. I. 
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Fig. 4 As in Fig. 3, but the scenario (b) in Fig. 1 is assumed. 
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3. The implications to the long baseline experiments 

In the present case with mass scale hierarchy the oscillation probability 
P{y^ — > u e ) in vacuum is given by 

P(u, - v e ) = 4\U e3 \ 2 \U, 3 \ 2 sin 2 = 4 sin 2 2^ sin 2 , 

so that we observe that the factor 4|f/ e 3| 2 |f/ M 3| 2 corresponds to sin 2 29 in the two 
flavor framework. For each m 2 in Figs. 3 and 4 the maximal allowed values of 
4|[y e3 | 2 |[/ M3 | 2 have been evaluated and are shown in Figs. 5(a) and 5(b) in the two 
flavor plot (sin 2 29, Am 2 ), respectively, where the regions which can be probed by 
K2K 23 , MINOS 21 and KamLAND 26 ' 25 are also given. The only promising channel 
for the K2K experiment is <-> disappearance, while MINOS might have a 
chance to see — > u e appearance and KamLAND may be also able to observe 
V e <-> u e disappearance. 
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Fig. 5 

Shadowed regions are allowed from the Superkamiokande atmospheric 
neutrino data for the scenario (a) (Fig. 5(a)) or (b) (Fig. 5(b)) in Fig. 
1, respectively. The region which can be probed by the long baseline 
experiments (K2K, MINOS, KamLAND) is the right side of each line. 
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4. The threefold maximal mixing model 221413 

Using (|1|) and (|3]), and redefining v e — ► e J<5 z/ e we have 
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where 6 m is the mixing angle in matter given by 
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as in the two flavor case 20,28 , and 



B = V (AE31 cos 2# 13 - A) 2 + (A£ 3 i sin 2fl 13 ) 2 . 

If we assume that the density of the Earth is approximately constant, then 
can be integrated as 
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The oscillation probability is given by 
P{v e — > u e ) = 1 — sin 2 29m sin 2 — - — 
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One of the interesting applications is the behaviors of the oscillation probability 
(||) in the three flavor maximal mixing model 22 ' 14 ' 13 , which is obtained by putting 
#12 = #23 = tt/4, #13 = sin" 1 (l/y/3j, 5 = tt/2 in the present parametrization. If 
|A£?3i| <C A then we see from (|5]) that the mixing angle 9m becomes tt/2, so we 
get from (|j) 

P{y e ^v e ) ~ 1 

P(z/„->z/ M ) ~ 1 -sin 2 2# 23 sin 2 | = 1 -sin 2 |, (6) 

where we have used the fact 

$ ~ (1 + cos2# 13 ) = c 2 l3 AE 31 L = -AE 31 L. 

(0) indicates that the mixing becomes pure «-> i/ T with maximal mixing and 
this is why a fit of this model to the data becomes reasonably good for lower 
Am^ 9 . On the other hand, if |A_E 31 | ^> A, then 6 M becomes 9\ 3 and we find 



again from 



P{v e ~^v e ) ~ 1 - ~ sin 2 29 l3 = - 



P( //,., — //,, ) ~ I - Lsi 3 sin 2 29 13 - sin 2 26 23 (s 2 3 sin 2 + ^c? 3 



5 1 2 AL 
sin . 

9 3 3 



where we have used 



AAL 2 
$ ~ (1 + cos 29 13 ) = c\ 3 AL = -AL. 

Possibility with larger Am 31 is excluded by the CHOOZ result 2 , but a fit of this 
model to the atmospheric neutrino data themselves is reasonably good because of 
the factor sin 2 (AL/3) accounts for the zenith angle dependence to some extent 6,9 . 
These two behaviors can be confirmed numerically as is shown in Fig. 6 where the 
value of Xchooz i s & l so plotted. Thus we see the analytic formulae for constant 
density explain the qualitative behaviors of the oscillation probability well. The 
threefold maximal mixing model fits to the data of atmospheric neutrinos and the 
CHOOZ experiment for m 2 ~8x 10~ 4 eV 2 9 . 
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Fig. 6 

X 2 fit as a function of Am 2 to 
the Superkamiokande atmospheric 
data for the threefold maxi- 
mal mixing model. The solid 
line includes both R and up- 
down asymmetries whereas the 
dashed line includes only the 
up-down asymmetries (# of de- 
grees of freedom = 6-1=5 for 
R+Y and =4-1=3 for Y only). 
The x 2 for the CHOOZ reac- 
1 tor data (dotted line; # of de- 
grees of freedom = 12-1=11) 
is also shown. 

5. Conclusions 

I have presented a result on the three flavor analysis of the Superkamiokande 
atmospheric neutrino data for 414 days. If the mass squared difference is smaller 
than 1 x 10 _3 eV 2 then there can be large mixing between and v e . It was shown 
that K2K will be able to see only <-> disappearance, while there may be a 
chance for MINOS (KamLAND) to see — > v e appearance {p e <-> v e disappear- 
ance), respectively. The threefold maximal mixing model gives a reasonable fit to 
the atmospheric neutrino and CHOOZ data for Am 2 ~ lx lCT 3 eV 2 . It should be 
emphasized that the matter effects are important in the analysis of atmospheric 
neutrinos. 

Acknowledgement 

The author would like to thank H. Minakata and E. Lisi for discussions, 
R. Foot and R.R. Volkas for collaboration and discussions. This research was 
supported in part by a Grant-in-Aid for Scientific Research of the Ministry of 
Education, Science and Culture, #09045036, #10140221, #10640280. 



CM 



40 
35 
30 
25 
20 
15 
10 
5 





0.0001 





! %itm(R+Y)_ 




xitm(Y) 




X 2 _ 

■^chooz 









0.001 0.01 

Am 2 / eV 2 



12 



6. References 



1. B. Ackar et al, 1995, Nucl. Phys. B434, 503. 

2. M. Apollonio et al, 1998, Phys. Lett. B420, 397. 

3. See, e.g., J.N. Bahcall, R. Davis, Jr., P. Parker, A. Smirnov, R. Ulrich eds., 
1994, SOLAR NEUTRINOS: the first thirty years Reading, Mass., Addison- 
Wesley and references therein. 

4. V. Barger, K. Whisnant, S. Pakvasa and R.J.N. Phillips, 1980, Phys. Rev. 
D22, 2718. 

5. G.L. Fogli, E. Lisi, D. Montanino and G. Scioscia, 1997, Phys. Rev. D55, 
4385. 

6. G. L. Fogli, E. Lisi, A. Marrone, and D. Montanino, 1998, Phys. Lett. B425, 
341. 

7. G. L. Fogli, E. Lisi, A. Marrone, and G. Scioscia, 1998, [hep-ph/ 9808205 . 



8. R. Foot, R. R. Volkas and O. Yasuda, 1998, Phys. Rev. D58, 13006. 

9. R. Foot, R. R. Volkas and O. Yasuda, 1998, Phys. Lett. B433, 82. 

10. Y. Fukuda et al, 1994, Phys. Lett. B335, 237. 

11. Y. Fukuda et al, 1998, Phys. Lett. B433, 9; 1998, [hep-ex/9805006 . 

12. Y. Fukuda et al, 1998, [hep-ex/9807003 . 

13. C. Giunti, C. W. Kim and J. D. Kim, 1995, Phys. Lett. B352, 357. 

14. P. F. Harrison, D. H. Perkins and W. G. Scott, 1995, Phys. Lett. B349, 137; 
1997, Phys. Lett. B396, 186. 

15. K. S. Hirata et al, 1992, Phys. Lett. B280, 146. 

16. ICARUS experiment, http : / / www. aq uila. inf n. it / icar us / . 



17. T. Kajita, Talk at XVIII International Conference on Neutrino Physics and 
Astrophysics, June, 1998, Takayama, Japan 
( [http: / / www-sk.icrr.u-tokyo.ac.jp/ nu98 / scan /063/ 



18. E. Lisi, Talk at New Era in Neutrino Physics, June, 1998, Tokyo Metropolitan 
University, Japan 

( |http://musashi. phys. metro-u.ac.jp/neutrino/junel2/e-lisi/ 



19. Z. Maki, M. Nakagawa and S. Sakata, 1962, Prog. Theor. Phys. 28, 870. 

20. S. P. Mikheyev and A. Smirnov, 1986, Nuovo Cim. 9C, 17. 

21. MINOS experiment, http: / / www.hcp.anl.gov / N DK / HypcrTcxt / numi.html . 

22. R. N. Mohapatra and S. Nussinov, 1995, Phys. Lett. B346, 75. 

23. K. Nishikawa, preprint INS-Rep.924 (1992); Talk at XVIII International Con- 
ference on Neutrino Physics and Astrophysics, June, 1998, Takayama, Japan 
( [http: / / www-sk.icrr.u-tokyo.ac.jp/ nu98 / scan /093 /| ) . 

24. Particle Data Group, 1998, Eur. Phys. J. C3, 1. 



13 



25. F. Suekane, Talk at New Era in Neutrino Physics, June, 1998, Tokyo Metropoli- 
tan University, Japan 

( |http://musashi.phys.metro-u.ac.jp/neutrino/junel2/f-suekane/| ). 

26. A. Suzuki, Talk at XVIII International Conference on Neutrino Physics and 
Astrophysics, June, 1998, Takayama, Japan 

( tittp : / / www-sk.icrr . u-tokyo .ac.jp/ nu98 / scan/ (J83/| ) . 

27. G. S. Vidyakin et al, 1994, JETP Lett. 59, 390. 

28. L. Wolfenstein, 1978, Phys. Rev. D17, 2369. 



